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Exact solution of qubit decoherence models by a transfer matrix method 
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We present a new method for the solution of the behavior of an ensemble of qubits in a random 
time- dependent external field. The forward evolution in time is governed by a transfer matrix. The 
elements of this matrix determine the various decoherence times. The method provides an exact 
solution in cases where the noise is piecewise constant in time. We show that it applies, for example, 
to a realistic model of decoherence of electron spins in semiconductors. Results are obtained for the 
non-perturbative regimes of the models, and we see a transition from weak relaxation to overdamped 
behavior as a function of noise anisotropy. 

PACS numbers: 68.65.Fg, 85. 35. Be, 03.67.Lx, 76.30.Pk 



I. INTRODUCTION 

In the last decade, advances in fabrication and exper- 
imentation have made it possible to observe and control 
quantum effects in individual systems that could previ- 
ously be observed only in the aggregate. One important 
condition for these advances has been to minimize, or 
at least to manage, decoherence. Quantum systems are 
inevitably subject to random influences that destroy co- 
herence, and classical behavior takes over. Only if this 
can be avoided can we hope to advance the technology 
of quantum information processing, and theoretical un- 
derstanding can help to achieve this. Theoretical calcu- 
lations generally fall into three categories: eactly soluble 
examples, approximation methods, and numerical com- 
putations. The theory of quantum decoherence has been 
notable for its richness in the last two of these three, but 
exactly soluble models have been few and far between, 
though some exceptions may be noted [1],[2],[3]. This 
paper goes some way towards changing this situation. 

Most treatments of decoherence begin with two cou- 
pled quantum objects: the system to be observed, and 
the unmeasured environment to which the system is cou- 
pled. The system eventually loses its coherence due 
to that coupling. The most natural way to treat the 
problem is to formulate a master equation for the re- 
duced density matrix, which incorporates the system- 
environment coupling, and also some aspects of the en- 
vironment's dynamics. This approach can in principle 
treat the "back-action" of the system on the environ- 
ment, though this is often difhcult in practice. The 
classic example is the Caldeira-Leggett model [4] 

An alternative formulation is to treat the environment 
as a fixed source of random noise. In this approach it is 
not possible to take into account back-action. However, 
there are many physical situations in which back-action 
is not important. Treating the environment as fixed but 
random is then usually more convenient. This approach 
will be used in this paper. 

In quantum computation, a large array of qubits - two- 
level quantum systems - needs to maintain coherence. A 
necessary requirement for this is to maintain quantum 
coherence of a single qubit, and most discussions begin 



at this level. The restriction to a two-state quantum 
system allows simplifications in the formalism. The pur- 
pose of this paper is to show that the decoherence of a 
two-state system subject to random noise can, under cer- 
tain conditions, be treated by a powerful transfer matrix 
method analogous to that of statistical mechanics, allow- 
ing us to solve a number of cases exactly. Furthermore, 
these solutions apply to the non-perturbative situation in 
which the coupling to the environment is not smaller than 
the level separation of the qubit. This difhcult regime 
is of great importance since some quantum computation 
schemes, most notably holonomic quantum computing, 
propose to operate in this parameter range [5] . 

In the next section we introduce and solve a class of 
white-noise-type models that can be solved by a relatively 
straightforward application of the method. We discuss a 
simple example, and compare the exact results to pertur- 
bation theory. The new results in the non-perturbative 
regime are also explored. In the third section we extend 
the method to Markovian models that can also be solved, 
this time by a slightly more intricate form of the trans- 
fer matrix. In the final section we discuss the possible 
extensions and ultimate limitations of the method. 



II. WHITE-NOISE-LIKE MODELS 



A. Definition 

Our two-level system will be described by the Hamil- 
tonian 



H = -Boa, - b{t) ■ a, 



(1) 



so that the energy separation is 2Bq and b{t) is a random 
function. The system is prepared at time t = with a 
2x2 density matrix p{t = 0) = pq. The measurable 
quantities are 



M{t^tf)^ / Vm]Vm] Tr [poa{t)] 



(2) 



Here axjCy, and a, are the Pauli matrices, which to- 
gether with the unit matrix represent a complete set 



2 



of observables in the qubit space. We work in the 
Heisenberg representation, so that a{t) satisfies i da/dt = 
[a,H]. J 'P\b{t)]T>[b{t)] indicates a sum over all possible 
functions b{t) on the interval t G [0,tf], with the appro- 
priate probability T'[b{t)]. The specification of V[b{t)] 
defines a model in the general category of 'qubits subject 
to random time-dependent forces. 

Our notation is of course motivated by spin qubits. 
However, it should be obvious that the method and all 
results apply to any two-level system. In particular, Bq 
and b need not be magnetic fields. 

We now consider the following form of P[b{t)]. Let b{t) 
be piecewise constant with discontinuous jumps at regu- 
lar intervals of length t. Furthermore, let b{t) be com- 
pletely independent between time intervals, i.e., b{t) = bi 
and B{t) = Boz + h = Bi for < i < r, b{t) = $2 and 
B{t) = Boz -I- &2 for T < i < 2t, etc. , with no corre- 
lation between the difi^erent b^. Finally, each 6j has the 
same probability distribution Let there be m time 

intervals so that the final time when a measurement is 
performed is = mr. For this model, then, we have 



jv[b{tmm]F^{t) 
= j d%p{h) ■ ■ ■ 



d'bmP{bm)F b{t) 



F b{t) 



for the average of any functional F b{t) . P{b) > and 

J d^b P{b) = 1. Any model with the piecewise constant b 
on equally-spaced intervals and the product form for the 
probability is exactly soluble, as we now show. We shall 
also give a physical example below, in order to indicate 
that the model is not artificial. 



Defining 



B. General solution 



Hi = —Bi ■ a = —BoCTz — bi - a 



and 



Ui = e-'"^\ 



the solution for a{tf) is 



aitf) = a{mT) = U+ ■ ■ ■ U+U+a{0)UiU2 ■■■Um- (3) 



Once we have a{tf) we can apply the initial conditions 
and obtain results for any observable according to Eq. 
2. Because of the statistical independence of the values 
of b at diflferent intervals, the averaging can be done for 
each correlation time individually, so aij transforms in 
the time r into 



aij{r) = I d%P{h){U+).^aki{0) {Ur\ 



[exp {iHiT%^ 



exp 



.<?fe;(0) 



Ecos{Bit) — iBi ■ asm{BiT) 



Ecos{Bit) + iBi ■ CTsin(Bir) ^^((O), 

J ij 



r 



and summation over repeated indices is implied. The 
identity matrix E always transforms into itself: E{t) = 
E. We have used the identity ex.p{iX ■ at) = cos Xt + iX ■ 



asmXt, valid for any vector X. X = X/ 



X 



x/x. 



The key to the method is to compute and then iterate 
the superoperator | [C^i(T)]y. Performing the 

average and writing the results in matrix form, we find 

= 7o?(0)-F^7yC7ia(0)aj+i^/i [<T(0)(Ti - 0-^^(0)] , 



where the integrals lo, li and 7^ are defined as 
7o = y P{b) cos^ {BT)d^b 

Ii = j P{b)Bi sin(Br) cos (Br) d^b 



(4) 



(5) 



P{b)BiBj sm"^ {BT)dH. 



Here B = Bq + b, B^ ^ b^/B, By = by/B, and B^_ = 

{bz + Bq) /B. (Note that Bq may be defined so that Ix = 
ly = 0. In order to exhibit the symmetry of the following 
expressions, however, we shall retain these quantitities.) 
Using the properties of the Pauli matrices to simplify the 
expressions, we find the three operators: 
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CTx(t) = (/O + 4X - - IZZ) <Tx(0) + 2 {I^y + I,) ay{Q) + 2 (7x, - ly) f7,(0) 

= {h - Ixx + lyy - hz) (Ty{Q) + 2 {hy - h) CTx(O) + 2 {ly^ + 4) a,(0) (6) 

fT.(T) = (7o - 7^^ - 7yy + 7,,) (7,(0) + 2 (7^, + ly) (7^(0) + 2 (7^,, - 7^) ay{Q). 



We may write this in a compact fashion by defining the 
transfer matrix T by 



o"x(r)\ 


1 ^xx 


T 

J-xy 


Txz\ 


fcTxiO) 






T 

^yy 


Ty. 


ayiO) 


'yz{r)) 




T 




Vz{0) 



or 

a(r) = T a(0), 

for short - it is understood that in expressions of this kind 
CT is a column vector. The detailed expressions for the 
matrix elements are: 







7o + 


^xx ^yy 




T 




/o- 


^xx ~i~ ^yy 


-^zz 


T 




/o- 


^xx ^yy 


+ Izz 


T 

J-xy 




^■^xy 


+ 2Iz 




T 

J-yx 




IT 


-2Iz 




T 

xz 




^^xz 


-2Iy 




T 

zx 




27a;2 


+ 2Iy 




T 

J-yz 




IT 


+ 27, 




T 

^zy 




IT 


-27,. 





(7) 



T is an example of an exactly computable quantum dy- 
namical map. For a definition and general discussion of 
such maps, see Ref. [6]. At tf, we find 



a{tf) = a(mT) 



j dH2P{b2)x 



(8) 



SO this process can be iterated, and we have 

&{tf)= &{mT) = ct(0). 

T is real but generally not symmetric. It oper- 
ates in the complex three-dimensional vector space of 
superoperators. Each vector {ax,ay,az) in the space is 
associated with the operator a ■ a. 
Diagonalizing T, we find 

D = RTR-^ 

where D = diag(di, d2, d^) and the rows of 7? are the 
eigenvectors of T. Since T is real, the eigenvalues are 



the roots of a cubic equation with real coefficients, so at 

least one of the di is real and the other two must either 
be real, or they must be complex conjugates. Iterating, 
we have 

a{tf) = ainiT) 
= a{0) 
= {R-^DR)"'a{0) 
= R-^D'^R a(0), 

with 

\dir 

T)™ = I {d2T 

{dzY 

Recalling that m = tf/r and using index notation, the 
solution becomes 



(9) 



This equation, together with Eqs. 5, 7, and 2 constitute 
an exact solution for all P{b). The magnitudes of the 
eigenvalues dj determine the relaxation times through 
the relation 



1 



1 



■ln\dj\, 



(10) 



while the transformation matrices 7? contain the informa- 
tion that allows an experimenter to prepare and measure 
a particular relaxation time. 



C. Realistic Example 

A relatively simple example of current interest con- 
cerns the spin dynamics of an electron in a two- 
dimensional system such as an appropriately constructed 
semiconductor heterostructure. Let the system be sub- 
ject to an electric field E — Ez and an applied uniform 
time-independent magnetic field BqZ perpendicular to 
the conducting layer. Due to spin-orbit couping, there 
is also an effective magnetic field acting on the spin, ac- 
cording to 6 = ap X E, where the electron momentum p 
lies in the x — y plane, a is a constant and \'p\ = Pf, the 
Fermi momentum. The electron scatters, causing p and 
therefore 6 to be a random function of time. The mag- 
nitudes \p\ and therefore b = bo are fixed, however. If 
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we take the scattering events to occur at equally-spaced 
intervals of length r, now understood as a time related 
to the momentum scattering time Tp, then this situation 
falls into the class of models considered in this paper. 
The assumption of equally-spaced intervals is of course 
an approximation - the actual length of the intervals is 
random, and is governed by a Poisson distribution. The 
effective magnetic field is called the Rashba field [7] , and 
the spin relaxation mechanism is known as the DP or 
D'yakonov-Perel' mechanism [8]. 

Hence the Hamiltonian for the spin degree of freedom 
is precisely of the form shown in Eq.l, and in this case b 
lies in the x — y plane. The average over h is performed 
by settting h — 6o(cos 0, sin i^, 0) and averaging uniformly 
over (j). B = Egz + b. In the physical situation of interest, 
the environmental correlation time r is usually much less 
than the spin relaxation times and also much less than 
1/Bq. Now set Ix = ly = 0, and note that symmetry 
in the x - y plane also gives Ixy = lyx = hz = hx = 
lyz = Izy = 0. For this particular case I^^ = since 
6^ = 0. but it is interesting to keep this quantity for the 
time being. The transfer matrix simplifies considerably: 



perturbative case first: 



^xx — -^0 ~t~ ^xx ^yy 
^yy — '^O ^xx ~1~ ^yy 
Tzz — -^0 ^xx ^yy 



(11) 



'-xy 



•-yx 



2L 



T = T = T = T =fl 

-^xz — ^zx — ^vz — ^zv — 



The eigenvalues of T are dz = Tzz = Iq — Ixx — lyy + 

Izz, and dx,y =h- Izz ± t^47f - {I^x - lyyf ■ In the 

Appendix we show that these eigenvalues satisfy Ids] < 1 
and \d±\ < 1, which guarantees that all solutions are 
exponentially decaying. 

The connection to dccoherence times is obtained by 
considering the limit where tf >> I/Bq >> r and tf >> 

1/ b » T. Expanding the integrals for small r, we find 

« 1 - 2t2 I P(6) {bl + bl) d'b=l- 2t2 (1+ 6|) 

(12) 

dx,y ~1 - / ^(^) (-^0 + fe^) d''b+ 

-T^ J P{b)Bl dH ± 2isjT^Bl-{Ixx-Iyyf 1^ 

=l-2T'B^,-r'(bl+¥y + 2M) + 

±2i^T^Bi-^^r^(¥x--^y (13) 

The result for dz is independent of the ratio b /Bq, but 
the dx,y are sensitive to it. Most experiments are done 
in the regime where b « Bq, and we consider this 



1 - 2t^B^ - (bl + bl + 2bl^ ± 2iTBo + 0{t^ 
l-r2(6| + 6|+26f) (14) 



If the system is prepared in an energy eigenstate, (in spin 
language, along the z-direction), then it remains along 
that direction since Txz = Tyz = 0, and we find 



(a,(t/)) = (d3)*^/^cT,(0) 



1-2t^ 



{bl + bl)]''^\az{0)). (15) 



This is finite in the limit of small 1 /b and r if the combi- 
nation T (bx + by^ approaches a finite limit, and we have 

{az{tf))^exp{-tf/T,){aM), (16) 
with the energy relaxation time Ti given by 

l-=2T(¥x+^)=2blT. (17) 

If the system is prepared in an equal superposition of the 
energy eigenstates (in spin language, in the x-y plane), 
then there arc oscillations of and (jy owing to the fact 
that the d^^y are not real, at an angular frequency given 
by <P/t, where (f) = |tan~^ (Imrf+/Red+)| . There is also 
exponential decay given by 

\{c^x,y{tf))\^^^Pi-tf/T2)\{cTx,ym\ (18) 

with the phase relaxation time T2 given by 

^=r(6|+62+26f)=6^r. (19) 

Note that these results satisfy_t;he constraint that in 
the isotropic {Bq 0) limit b^. = by = b^ we must 
have Ti = T2, since the distinction between longitudi- 
nal and transverse relaxation ceases to have any mean- 
ing. Incidentally, the perturbation theory results for the 
same model except that the time intervals follow Pois- 
son statistics are: 1/Ti = 4,Tp (bl + by^ and I/T2 = 

2Tp (bx + bl + 26^^ [9]. Thus the previous results can 
be translated over with the identification r = 2Tp. 

The rotation matrix R is block-diagonal for this sim- 
plified case: Rxz = Rzx = Ryz = Rzy = 0. This means 
that preparation and measurement of the spin in the 
2— direction will give Ti and preparation and measure- 
ment of the spin in any direction in the x — y plane will 
give T2. Preparation in any other direction will result 
in independent relaxations of the z— component and the 
magnitude of the projection of the spin onto the x — y 
plane. There is no cross-talk between the two relax- 
ations. 
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D. Comparison with perturbation theory 

The usual treatment of problems in this class is by per- 
turbation theory, often going under the name of Rodfield 
theory in this context. The theory has been extensively 
worked out because of its applications in nuclear mag- 
netic resonance. The condition for the validity of the 
perturbative approach is >> b{t) ^.t all times. A 
detailed treatment is given by Slichter [10]. The results 
for relaxation times are: 



J 1 



(20) 

(21) 



with Wo = 2So and 



kii{Lo) =2 ds e-''^'bi{t)b,{t + s). (22) 



To compute kii in the example of the previous subsection, 
we note that fc^x = kyy and k^z = and that we must 
average t uniformly over the interval (0, r) and average 
(f> uniformly over the interval (0, 27r) . For the piecewise 

constant b{t), we find 



bi{t)bi{t + s) =b^cos [(/)(t)] cos [(/)(t + s)] 



= <! bl/2, 
0, 



if 
if 
if 



s < -t 

-t<S<T-t 
S>T -t. 



(23) 



This result can be substituted into Eq. 22. In the limit 

that Bqt << 1, kxxi^^o) ~^ bgT, and similarly for kyy^LOo). 
Thus perturbation theory applied to our model yields 



l/Ti = 2&§T and I/T2 = l/2Ti, in agreement with Eqs. 
17 and 19. Although we have not treated them in de- 
tail, the longitudinal fluctuations described by the func- 
tion kzz that produce dephasing also have the correct 
perturbative limit [cf. Eqs. 19 and 20]. 



E. Non-perturbative regime 

In this subsection, we wish to investigate the overall 
behavior of the exact solution. There are four param- 
eters in the model: r, I/Bq, 1/ b , and r/. The regime 
tf < r is of no interest: the spin processes freely about 
the random field Bq +bi. So we focus on the regime 
tf » T, and the expression for the observables in Eq. 
9. We have already discussed the perturbative regime in 

which tf >> 1/ b >> 1/Bq >> t. Our solution is non- 
perturbative, and remains valid for arbitrary values of 



/Bo, 



T and Bqt. However, simple expressions are 



available only in relatively symmetric situations, so let us 
take b =bo and I^y = lyx = hz = hx = lyz = Izy = 0, 
but allow for the possibility that Ixx 7^ lyyi and also 
that bz ^ 0. We also assume the probability density 

P(6)posesses wcll-dcfincd moments up to the fourth or- 
der in 6, and that the odd moments vanish. We once 
again have dz = Tzz = /o - Ixx - lyy + hz- 

Since we have a closed-form solution, we can easily 
compute higher-order expansions. To see this, take 

r << 1 and Bqt << 1 in order to perform the trigono- 
metric integrals. Then to order r^, we have 



= --\n\dz\ = 2{bl + bl)T+\2(bl + bl) - - [2BI {bl + bl) + 2b% + 264 + 2blbl + blbl + blbi\ \ (24) 



The strength of the steady field enters the relaxation rate at order t^. 



For T2 there are two possibilities. For 1^ > \Ixx — Iyy \ j 

All - {Ixx- lyy f, \dxf = {lo ' hzf + 4/f - [hx ' lyyf and 



dx has an imaginary part, dx,y = Iq — Izz ± 



= --In 141 = --ln\dy\ 

±2 T T 



(25) 



= r{bl + bl + 2bi) + (bl + bl + 2bl) ' Isgfe^ + i (5^ + 6^ + 2bi + 2blbl + ^b\bl) - '-{^ 



1 /_ _x2 



bl - bl 



again good to order r**. This corresponds to a damped 
oscillatory solution. For 1^ < \Ixx — Iyy \ 1 the dx,y are 
real and the solution is overdamped. All three eigenval- 
ues of T are real, and so there arc three relaxation times. 
Ti remains as above, and, to order r^, we find 



^ = (62 + 62 + 262) r ± y(62-62)'-4BoV. 
Here the spin relaxes so fast that it does not have time 
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to precess about the external field. 



(a) 



Of particular interest is the fact that this abrupt tran- 
sition from oscillating to overdamped behavior is driven 
by the anisotropy of the noise. Physically, strong highly 
anisotropic noise quickly erases the distinction between 
the clockwise and counterclowise directions about the z- 
axis, and the spin does not know which way to precess. 



The most interesting issue that is accessible for the 
exact solution is the crossover from the weak-noise regime 

( b « Bq) to the strong-noise ( b » Bq) regime. 



We first loqk_ at the case of purely transverse noise: 

6^ = 6^ and 5^ — 0. The probability distribution is 
again the uniform one for the angle 9 = ta,\i~^{by/bx) 
with b"^ + by = 6q, a constant. In Fig. 1 we plot the 

normalized rates l/Tib'^r and l/T2b'^T as a function of 
/Bq for three values of Bqt : Bqt = 0.1, Bqt = 1 
and Bqt = 10. First note that the correct perturbative 
limit is obtained when b /Bq < 1 and Bqt < 1 in Fig. 
1(a). The standard transverse-noise relation 1/Ti = 
2/T2 holds over a fairly wide range: roughly 6 /Bq < 3, 
after which 1/Ti > 2/T2. Stronger noise causes 1/Ti 
to be a super-quadratic function of 6^, while the usual 
physical picture of dephasing (I/T2) as due to a random 
walk in the phase variable appears to hold for all values 

of b /Bq. As Bqt increases, Fig. 1(b) and 1(c), we find 
that the relaxation rates decrease. This is not surprising: 
as T increases, k^x (^) = kyy(uj) become sharply-peaked 
functions of and the resonant processes responsible for 
the relaxation arc suppressed. As the rates decrease, the 
1/Ti = 2/12 holds over a braoder range of noise strength 

b /Bo. 



Secondly, we consider the case of isotropic noise: 6^ = 
by — b1 with a uniform distribution of the solid angle and 

6^ 



y 



= 5q, a constant. In Fig. 2 we plot the normal- 
ized rates l/Ti^^r and l/T2b'^T as a function of b /Bq 

for the same three values of B^t : Bqt — 0.1, Bqt = 1 
and Bqt = 10. First note that the correct perturbative 

limit is obtained when b /Bo < 1 and Bqt < 1 in Fig. 

2(a): the relation l/Ti = l/r2 holds quite well for all 

values of noise strength b /Bq, but both 1/Ti and l/r2 

are super-quadratic functions of b^. As Bqt increases, the 
dephasing rate 1 /T2 begins to dominate over the energy 
relaxation rate l/Ti. Figs. 2(b) and 2(c), the relaxation 
rates decrease for the same reasons as above. Inter- 
estingly, however, the l/Ti — I/T2 relation is violated, 
presumably again due to the varying relative weights of 
kxx{^) and kzz{io = 0) as r is increased. 
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FIG. 1: l/(ri,2FT) as functions of |b| /Bo for 6|=b| and 

b1 = Q with BoT = 0.1, 1 and 10, corresponding to (a), (b) 
and (c) respectively. Upper (red) lines are for 21 and lower 
(blue) lines for T2. 



III. CORRELATED-NOISE MODELS 



A. General solution 



The Hamiltonian remains the same: 



H = -Boaz - b{t) - a^-B-a, 
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FIG. 2: l/(ri,26l^T) with bj^ = as functions of |6| /Bo 

for bi=hy = hi with _Bot = 0.1, 1 and 10, corresponding to 
(a), (b) and (c) respectively. Lower (red) lines are for Ti and 
upper (blue) lines for T2. 



and the piecewise constant time dependence of h{t) is 
retained. Now, however, successive values of hi may 
be correlated. Thus we take the following model for 
P(6i, b2, ■■, bm) ■ the values of hi are independent except 
for neighboring intervals, so we have 



We shall make the natural assumption that P{b, h') is 
symmetric: P(6, 6') = P[b',b). 

Now (T,-, transforms in the first interval into 



Mr){b2) = j db,P{b,M) Mir)]^^ I [C/i(r)]y au 
= J dbiP{bi,b2)[expiHiT].j^ exp-iHiT 

= j dbiPibiM) 

X i?cos(i?iT) — iSi • (Tsin(_BiT 

X E coii{BiT) + iBi ■ S sa\{BiT) cfki, 
L J ij 

(26) 

where E is the identity matrix and summation over re- 
peated indices is implied. In matrix form, this is 



+ i ^ (&2) i^o-t - o-^ct) , 

i 

where the integrals Iq, E and Ej are defined as 
/o(6') J P{b,b') coa^ {Bt) db 
E{b')= / P{b,b')B,sm{BT) cos (Bt) db 



(27) 



/,y(&') = / P{b,b')B,Bj ain^ (Bt) db, 



analogously to the previous case. 

We may write the evolution equations in a com- 
pact fashion by defining /o(6) = cos^ (Bt), fi(b) = 
i?i sin(i3r) cos (i?T) , and fij{b') = BiBj am^{BT), and 

Txx{b) = /o(fo) + fxxib) - fyy{b) - fzz{b), ctc, as in Eq. 
7. Then we find 

a(62,T) = j dbiP{b2,bi)T{bi) cj{0), 

where a is to be interpreted as a column matrix. When 
this process is iterated we have 



P{bi,b2,..,b,n)^Tl, 



iP(6„6,+i) 



a(&3,2T) = / db2P{bMT{b2) <?(&2,t) 
db2P[bMT{b2) 
X j dhPihM) T{h) a(0), 
and finally 

a{mT) = [ dbm+i J dbjnPibm+i,bm)T{bm) • ■ • X 

db2P{hM)T(b2) I dhPihM) T{h) a{0). 
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This problem is exactly soluble when the function P{b, b') we have 
is separable: 



N 



(28) 



(Ti{mT) = j dbm+1 j dbrnPnAbrn+l)Pn,nibm)Tii^{brn) ■ 



j dbiPni 



{b2)PnAbi) Ti,j{b,) ajiO). 



This is of course not true for all functions P{b, b'). How- 
ever, if P{b, b') is continuous, it can be approximated uni- 
formly by such an expression [11]. Substituting and con- 
verting to index notation with repeated indices summed. This can now be written as a matrix product 



o-i(mr) = j 



dbm.+lPnm{bm+l)S{ 

"mi»)>("m-l>»m)'^("Tn-li*m))("n 



r 



\n2,h),{n„i2) J dblPn.ih) Ti,j{h) (7, (0), 

(291 



where we have defined the 3iV x SN matrix 

S(n',i),(n",j) = J db Pn'(b)Tij{b)pn"{b). (30) 

Diagonalizing S, we have the solution 

a,{mr)=Jdb'F^{b') (5-^) x 

J db F„-{b)T,k{b) afc(O) 

= I db'Fn{b'){R-s^ [diag{Sr-'] i?s 

X j dbF^,{h)Tjk{h) afc(O). (31) 

Here Rs is the matrix whose rows are the eigenvectors 
of 5, and diag{S) is the diagonal matrix with the eigen- 
values of 5* along the diagonal. The relaxation times are 
then obtained by the same arguments as in the previous 
section, except now they must be obtained by comput- 
ing the eigenvalues of S, not T. The matrices enclosing 
diag{S)"^~^ determine how to prepare and measure the 
system if one wishes to determine some particular relax- 
ation time. 



with < r < 1. We choose this model partly because 
of its simplicity. However, it also has some interesting 
physics. For r = it reduces to the previous example. 
For r = 1, there is suppression of backscattering at (j) — 
<p' = TT. This implies that, as r increases, the problem 
crosses over from a purely random walk in spin space to 
one in which the successive steps tend to be in the same 
direction. 



After performing the integrals over the 
find 



P{cl>, <f>') = ^[l + r cos {cP-cP')] 



and biz, we 



— (1 -h r cos(/icos0' -|- rsini^sin^') (33) 
2n 



so that the function is separable 



Pi<P,<t>') = Y^Pn{ct>)Pn{<l>') 



n=l 



(34) 



with 



B. Example 

In the first example, the scattering was purely s-wave: 
isotropic in </>. If the scattering can also occur in the 
p-wave chamiel. then the successive values of b are cor- 
related. An interesting example is a mixture of s-wave 
and p-wave given by 



PiW 



-j^, P2(0) = J-^COS(j}, psicj)) 



We now need to compute the 9x9 matrix 



'(m),( 



n'j) - j 



d(l) Pn{<P)Tij{(j))p„>{(j)). 



(35) 



(36) 



-» 1 r/T 7 \s:/. NT /I We are interested hcre in thc dccoherencc Hmit: b^r « 

^ ' ' " 2^0 ^ ~ °' ^ """^ ~ ->] ' (^^) BoT«l. The integrations over the angular functions 



9 



6 ^ 




0.2 0.4 0.6 0.8 1 



r 

FIG. 3: 1/(Ti,26^t) as functions of r. The upper (red) line is 
for Ti and the lower (blue) one for T2. 

are elementary. The nonzero elements of S are 

S(2.),(2x) = S(3y)^i3y) - (r/2) (l - 23^) - vhlr^ 

^(3.),(3x) = 5(2,),(2,) = (r/2) (1 - 2i?oV2) - 3r6gTV4 
'S'(iz),(iz) = 1 - 25or^ 

'S'(2z).(2z) = 'S'(32),(3z) = (''/2) (l - 26or^) 

= — '5'(li/),(l2;) = 2i?0''' 

'S'(22;),(2y) = 'S'(32:),(3i/) —S{2y),{1x) = ^'S'(3y),(3x) ^^0^ 

'5(22;), (3ij) = 5'(22;),(3ij) 'S'(2y),(3a:) = 5'(2y),(3a;) Th\T^ 

= S(iy)^Ciz) = 5(3^)^(12) = S'(l2),(32;) 

= 'S'(3z),(iy) = \/r&o5oT^/2 

'S'(l2;),(3z) = 'S'(3x),(l2) = -'S'(i.y)^(2z) = -S'(2y),(lz) 
= "'S'(lz),(3x) = -'S'(3z),(l2:) = '5'(lz),(2j,) 
= S(2z),{ly) = -V^fooT/2 

This matrix is easily diagonalized numerically, and the 
results are shown in Fig. 3. All eigenvalues are less 
than one in magnitude, as before. Six of these eigenval- 
ues are much less than 1, and correspond to transients. 
Two (one) give finite relaxation times when tf/r becomes 
large and correspond to T2 (Ti). 1/Ti increases from 26qT 
to about 5.9 6qT, and l/r2 increases from bpT to 2.5 h^r 
as X goes from to 1. This is to be expected, as decrease 
in the average angle of change of the field corresponds in- 
tuitively to a longer step in the random walk of the spin 
vector. 



IV. DISCUSSION 

The two basic features that make the basic white-noise- 
like model solvable is that (1) the noise Hamiltonian is 
piecewise constant with equal intervals and that (2) the 
system space is two-dimensional. These features limit 
the applicability of the model. 



Condition (1) is a fairly severe limitation: often, the 
effectiveness of the noise is related to fcii(ajo). In the basic 
model, this can only be varied by changing the parameter 
B{jT. The shape of kii{ujQ) is fixed in the white noise 
model, and can only be changed to a limited extent by 
going over to correlated- noise generalizations. Some very 
interesting problems, such as the effect of 1// noise on 
qubits [12], appear to be outside the scope of this method. 
It may be possible to allow for a distribution of intervals 
without sacrificing solvability - the multiple time integral 
thus obtained, though complicated, may be susceptible 
to integral transform techniques. 

Condition (2) is not easily jettisoned. The three Pauli 
matrices generate the 3-dimensional algebra of SU{2), 
T generates orbits in this algebra, and this implies that 
T is a 3 X 3 matrix. The transformations of a three- 
dimensional Hilbert space (for a qtrit) belong to SU{3), 
and the algebra is an 8-dimensional space. Hence the 
transfer T matrix for this problem would be 8 x 8. 
There may be special situations in which there is an 
exact solution. The solution of the two-dimensional 
Ising model depends on diagonalizing a 2^ x 2^ transfer 
matrix, where N is the number of spins in the trans- 
verse direction. The corresponding generalization of the 
present problem would require the diagonalization of an 
N"^ — 1 X — 1 matrix for an iV-level system. — 1 
is the number of generators of SU{N). Such a general- 
ization is likely to require some new ideas. 

Despite these limitations, the method is new and sus- 
ceptible to development in different directions. It should 
serve as an important tool in the study of decoherence in 
the future. 

We would like to acknowledge useful converstaions 
with C. Tahan and M. Friesen, and the support of the 
NSF ITR program. 



V. APPENDIX 



Here we show that the eigenvalues satisfy jdsj < 1 and 
\d±\ < 1. 

Theorem: |d±| < 1. 
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Proof: 



Jo - ± i V - - lyy) 



= I / P(6) [cos2(Br) - B,sm^{BT) 



1 2 



P{b)B^ sm{BT) cos {Bt) d^b 



< j cos^ {Bt) - 



B^sm^{BT) P{b)dH+ 



B^ sin(Sr) cos {Bt) P{b)d^b. 



where the final inequahty is an application of the Schwarz 



inequality. Since 



B, 



< 1, we also have 



|d±r</ cos^{Bt)- B^sm.{BT) P{b) d^b+ 

2 



+ 4 



/ 



B^ sin(Sr) cos (Br) P{b) dH 

cos^ {Bt) + B^ sin2(Sr)j ^ P{b) d^b 

< j [cog'{BT) +sm^{BT)Y P{b) d'^b 

= j P{b) d^b 
=1. 



Theorem: l^sl < 1. 
Proof: 



ds /q ^xx ^yy ~t~ ^zz 

P{b) [cos' {Bt) + -Bl~ Bl) siu^Bt) 



d% 



d% 



< j P{b) cos2(Bt) +Bf sin2(B 

< j P{b) [cos' {Bt)+ sin' (Bt)] d^b 



and recaUing that B' + B' + By = 1, we also have 



d3 = J P{b) [cos2(Br) + -Bl- sin2(Br)] 
= J P{b) {cos2(St) + [1 - 2 ^ 52^] sin2(Br)} d 
= 1 - P{b) cos' {Bt) + sin' {Bt) (b^ + B^^j d^b 
> 1 - 2 P{b) [cos2(St) + sin2(BT)] d^b 



Hence Ids] < 1. 
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